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ABSTRACT_ 

In a context of ever-growing electricity consumption and need for less polluting 
sources of energy, salinity gradient power (SGP) based on osmosis is a promising 
technology. Salinity difference between two solutions separated by a semi-permeable 
membrane leads to the pressure increase. The aim of this study is to find the critical 
permeability threshold of a membrane for the dimensioning an osmotic power plant. 
Using Spiegler-Kedem equations, the various fluxes across the membrane have been 
calculated, and delivered power is explicitly derived in terms of system parameters. 
A necessary condition for economic viability is that its upper bound is larger than a 
critical threshold value below which osmotic power plant is not profitable. As it is 
directly proportional to membrane permeability, fixing the optimal membrane 
permeability value will in turn enable conceive more efficient membranes specifically 
made for osmotic energy production, as such membranes do not exist today. 
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1. Introduction 

The considerable and very fast growing energy 
consumption consecutive to the economic 
development of many countries is acutely raising 
the question of avoiding disastrous 
environmental consequences which could 
inevitably occur if only conventional fossil 
sources are used. To cope with such situation, a 
large effort has been oriented toward other 
sources such as less damaging renewable ones 
(Lewis, ef a/., 2011; Kumar, ef a/., 2011). Aside 
solar, wind, geothermal and hydraulic sources, 
they also include other less evident ones which 
may however represent interesting alternatives in 
specific situations. 

Such is the osmotic power which is the process of 
converting the pressure differential between 
water with high salinity and water with lower or 
no salinity into hydraulic pressure (Loeb, 1975; 


Mishra, 2013; Kho, 2010; The European 
Commission, 2004; Heifer, ef a/., 2013; Skilhagen, 
ef a/., 2012; Skilhagen, and Aaberg, 2012). The 
harnessing of this energy for conversion into 
power can be accomplished by means of 
Pressure Retarded Osmosis (PRO) (Kim, and 
Elimelech, 2013; Heifer, ef a/., 2014; Wang, ef a/., 
2012). This technique uses a semipermeable 
membrane to separate a less concentrated 
solution, or solvent, (for example, fresh water) 
from a more concentrated and pressurized 
solution (for example sea water), allowing the 
solvent to pass to the concentrated solution side 
(Post, 2009). The additional volume increases the 
pressure on this side, which can be depressurized 
Dy a hydro-turbine to produce power and 
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electricity (Kleiterp, 2012) .As seventy percent of 
Earth surface is covered with water, 97 percent 
of which is saltwater, the process created by 
mixing seawater with freshwater generates a 
resulting osmotic power which could serve as 
both renewable and consistent electricity 
source. While still in early stages, best estimates 
of the global production potential of osmotic 
power exceed 1,600 terawatt-hours, the 
equivalent of half of Europe entire energy 
demand. There are two primary sources for 
osmotic power: 1) natural occurrence where 
river water meets the sea water, and 2) merging 
two man-made water sources from processing 
plants 

Both methods can be viable but in 1), seawater 
averages 40 grams of salt/ liter + River Water, 
less power is provided than in than 2), where 
brine (from desalination) averages 60 grams of 
salt/ liter + treated water. The higher the salinity, 
the more free energy can be extracted, and the 
more power can be generated. Today osmotic 
power is a promising renewable energy source 
(RES) provided conversion factor from pressure 
differential can be made large enough (Dinger, 
et at., 2012; Brcein, et at., 2010; Straub, et at., 
2016). In the following, some elements on this 
question will be discussed in PRO case. In 
particular, the threshold value for permeability 
coefficient which characterizes membrane 
efficiency for viable economic application is 
determined. 

2. Osmosis Pressure Representation 

Potential osmotic pressure n, the maximum osmotic 
pressure in a solution separated from osmosed fluid by a 
selectively permeable membrane, is given by 

n — CiRT (1) 

with 7i the potential osmotic pressure (Pa), C the solute 
molar concentration (mol m“ 3 ), R = 8,314 J. mol -1 . K -1 , 
T solution temperature (°K), and I the particle number per 
entity. Given the solutions, D and F with respective 
osmotic pressures 7 Cd and jif, where D is draw solution 
(most concentrated), and F feed one (less concentrated), 
osmotic pressure difference An = 71 d - 71f between the two 
solutions is typically equal to 12 bars for water fluid, ie, 
12.10 5 Pa. Let AP = Pd - Pf with Pd and Pf the pressures of 
solutions D and F. With Spiegler-Kedem [18] model flux 
equations can be integrated and one gets volume and solute 
fluxes J v and J s flowing through the membrane with 
pressure retarded osmosis 

J v = A(aA7t - AP) ; J s = - (1 - E m )->(J V C d E m - CO (2) 


with J v , J s (m.s -1 ) the volume and solute fluxes across the 
membrane, cj the solute reflection coefficient, A(m.Pa _1 
s' 1 ) the membrane fluid permeability, and E m = exp - [(1 

- a) JvLm '] with co the solute permeability. Similarly in 
the support one gets 

J s = - (1 - Es) -1 (J v CiE s - Cf) (3) 

with Q, Ci and Cf the solute concentration in the draw 
solution, at the barrier-layer/support interface, and in the 
feed solution respectively, and E s = exp - [J v 8D e _1 ] with 8 
the support thickness. Equality of fluxes in the barrier layer 
and the support gives from (1,2,3) the non-dimensional 
equation 

a{(l-a)(l-E‘- < V + E p (l-E p ) -1 ) = (l-rJCl-EP)' 1 - (4) 

with a = 1- Cj/Cd, a the solute reflection coefficient, r = 
Cf/Cd, p = 8co/D e L, 8 the support thickness, L the barrier 
layer thickness, D e the effective diffusion coefficient of the 
solute, E = exp-X, X = F(aa - <AP>), F = ARTvCdL/©, v 
the salt stoichiometric coefficient and finally <AP> = (Pd 

- Pf)/RTvCd. For a given dimensionless hydrostatic- 
pressure difference <AP>, (4) is a transcendental equation 
in a (dimensionless concentration difference across the 
barrier layer) which determines the operating conditions of 
the osmotic plant for a given set of system parameters. 

3. Osmotic Power Plant Production 

The basic energy production system is 
composed of a compression unit which delivers 
pressurized salted water injected in a chamber 
with a filtering membrane across which a flux of 
salted water J s is crossing, see Figure 1. 
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Figure 1. Sketch of Osmotic Plant Balance 
Js is Salt Flux across the Membrane. 

The power produced per surface unit of installed 
membrane is given W = J m AP with W in Watts/m -2 , 
or else in non-dimensional form 

<W> = W/A(RTvC d ) 2 = <AP>(aa - <AP>) (5) 
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<W> is maximum when <AP> = aa/2 and is then 
equal to <W> = (aa) 2 /4. But this is not necessarily 
possible as a and <AP> are also linked by (4) which 
imposes a constraint on system coefficients. One 
effectively gets for z = aa the transcendental 
equation 

z{(l-a)( 1 -E(z) 1_(J ) _1 + E(z) p (l-E(z) p ) -1 } = a (l-rXl-ECz)*’)- 1 - a 2 (6) 

with E(z) = exp-.5Fz which relates “physical” z 
(concentration performance) to “technical” F (barrier 
quality). From this different limits can be evaluated 
for possible power output from the system (Yip, and 
Elimelech, 2012; Yip, and Elimelech, 2011; Lin et al ., 
2014). However simple ones are directly obtained 
from (6) for <W>. In the case of very large p for 
instance, (4) takes the very simple form 

Y = - (Ln(l-x))/x (7) 

where x = a[l-r((l-a) -1 ] and Y = a(l-a-r)F/2. In 
definition interval x e [0,1], Y is monotonically 
increasing from 1 to + oo, and normalized delivered 
power <W> = Kx 2 , with K = o 2 /4[l-r((l-o)-'l 2 , 
varies monotonically between 0 and K, showing that 
larger <W> corresponds to larger x. More generally 
parametric dependence of power output W can be 
obtained from (6) to get best parameter range. When 
returning to dimensional expression, W is upper 
bounded by W sup = KA(RTvCd) 2 , and a necessary 
condition for economic efficiency is that 

W sup = KA(RTvC d ) 2 > W crit (8) 

where W cr it is the threshold value above which the 
osmotic plant is viable. W sup is larger with larger draw 
concentration Cd, larger temperature T and larger A, 
which has to be determined to satisfy economic 
system efficiency demand, ie here A > A cr it = 
4W C ritA7i 2 . For usual figure of delivered power Wd = 
5Wm' 2 , and with typical Arc = 12.10 5 Pa for water, one 
gets, for instance, A > A cr i t = 1.38 mPa _1 s _1 . More 
generally the variation of A cr it vs An, sees Figure 2, 
indicates for Arc e [6.10 5 ,20.10 5 ] the “efficiency” 
interval A e [4.910 _12 ,5.410 _n ] which is still at 
technical limit today (Straub, et al., 2016; Zhang, and 
Chung, 2013). To cope with these economic operating 
constraints, extension of simple osmotic barrier effect 
have been recently considered (Chou, et al., 2013; 
Chou, et al., 2012; Yip, et al, 2010; Banchik, et al., 2014; 
Me Cutcheon, and Elimelech, 2007; Cath, et al., 2013; 
Chen, et al., 2016; She% et al., 2016 ; Chou, et al., 2013; 
Hickenbottom, etal., 2016; Dechadilok, and Deen, 2006). 


A coupled system with solar plant will be discussed 
elsewhere. 


[— Values of A depending on the osmotic pressure 
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Figure 2 . Permeability Coefficient vs Osmotic 
Pressure Difference. 

4. Conclusion 

From analysis of equations representing osmotic 
physical phenomenon, it is shown that in Pressure 
Retarded Osmosis case, the harnessing of 
salinity-gradient energy taking place at the 
interface between waters of different salt 
concentration could provide an interesting and 
almost inexhaustible energy source if systems 
coefficients satisfy operating conditions which 
have been explicitly stated within Spiegler- 
Kedem model. However, even in optimum mode 
operating case, satisfaction of economic 
viability condition expressed by fixed specific 
membrane power output Wthres = 5 Wm n2 is not 
always met. Aside theoretical research on 
optimizing system operating mode, this 
weakness is urging further study of membrane 
physical properties, such as porosity and 
tortuosity pore length, to design most efficient 
hydraulic permeability of the barrier layer. 
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